Force generation in small ensembles of Brownian motors 
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The motihty of certain gram-negative bacteria is mediated by retraction of type IV pili surface 
filaments, which are essential for infectivity. The retraction is powered by a strong molecular motor 
protein, PilT, producing very high forces that can exceed 150 pN. The molecular details of the 
motor mechanism are still largely unknown, while other features have been identified, such as the 
ring-shaped protein structure of the PilT motor. The surprisingly high forces generated by the 
PilT system motivate a model investigation of the generation of large forces in molecular motors. 
We propose a simple model, involving a small ensemble of motor subunits interacting through the 
deformations on a circular backbone with finite stiffness. The model describes the motor subunits 
in terms of diffusing particles in an asymmetric, time-dependent binding potential (flashing ratchet 
potential), roughly corresponding to the ATP hydrolysis cycle. We compute force- velocity relations 
in a subset of the parameter space and explore how the maximum force (stall force) is determined by 
stiffness, binding strength, ensemble size, and degree of asymmetry. We identify two qualitatively 
different regimes of operation depending on the relation between ensemble size and asymmetry. 
In the transition between these two regimes, the stall force depends nonlinearly on the number of 
motor subunits. Compared to its constituents without interactions, we find higher efficiency and 
qualitatively different force- velocity relations. The model captures several of the qualitative features 
obtained in experiments on pilus retraction forces, such as roughly constant velocity at low applied 
forces and insensitivity in the stall force to changes in the ATP concentration. 

PACS numbers: 87.16.-b,05.40.-a 



I. INTRODUCTION 

Recent experimental progress has enabled remarkable 
quantitative measurement of biological processes on the 
single-molecule level 0] . One example is the biomechan- 
ics of force generation by molecular machines such as 
kinesin, myosin, and dynein This has stimulated 

considerable modeling activity in order to analyze the ex- 
periments H, IH lOl ' In this paper we are inspired by 
another motor protein, called PilT Q , which has interest- 
ing properties; e.g., it is the strongest known molecular 
motor 

The PilT motor is responsible for the retraction of cer- 
tain bacterial surface filaments, and the velocities and 
forces generated during retraction have been measured 
in a scries of laser tweezers experiments |^ |^ ^3 ■ Theo- 
retical analysis of the retraction data has revealed inter- 
esting information about the underlying retraction mech- 
anism . It is also of interest to study the question of 
how large the generated forces can be, given the energy 
and length scales relevant to the PilT motor, and what 
features are important for generation of large forces. We 
will address this question through a simple ratchet model. 
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which is inspired by known experimental facts of the PilT 
system. 

Ratchet models of particles in fluctuating potentials 
are commonly used in theoretical studies of molecular 
motors Single-particle models have been used to 

study kinetics of ATP consumption in molecular motors 
[T^ and to describe the kinetics of kinesin ^3 ■ Models of 
particles in ratchet potentials have also been employed to 
describe collective effects in larg e ensembles of interacting 
motors B El 111 El El El US- Finite ensembles of 
Brownian particles have been studied to some extent in 
the context of two-headed motor proteins 12^123 . 0, 
I25I I and to describe the bacterial flagella motor |26lj"^ 

Another approach, in which the motion of a molecular 
motor is described in terms of transitions between dis- 
crete chemical and conformational states, has also been 
generalized to the case of two interacting motor subunits 
jmAs is evident from Refs. [3 El El El El El El 
I21II22I I21 0, [2I I25 , interactions among several motors 
can lead to new and nontrivial behavior of the average 
velocity, which is not present if the interaction is turned 
off. It is natural to ask if this is the case for force gener- 
ation as well. 

In this paper, we investigate the behavior of a small 
ensemble of interacting processive Brownian motors and 
focus on the effect of the interaction on the generation of 
large forces. We aim at a prototypical, minimal model 
which captures certain features of the PilT system. The 
main input is the overall structures of the filament and 
PilT complex. When possible, we also use the exper- 
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imental situation to estimate model parameters, which 
should ideally be as few as possible. The model is a gen- 
eralization of the model of two elastically coupled mo- 
tors studied by Dan et al. to larger ensembles, but 
focuses on different properties. By varying the density 
of motor subunits and other parameters of the system, 
we explore the force production in different regimes. We 
compare with a single building block of our model to iden- 
tify the effect of the interactions and also compare with 
experimental results. Although the detailed connection 
between the model and the actual molecular retraction 
mechanism is speculative, the spirit of the model is best 
understood in light of the known facts about the PilT sys- 
tem. Therefore, we will briefly review some facts about 
pilus retraction before introducing the model. 

Type IV pili are surface filaments crucial for the ini- 
tial adherence of certain gram-negative bacteria to tar- 
get host cells, DNA uptake, cell signaling, and bacte- 
rial motility Each filament consists of thousands of 
pilin monomers that polymerize to a helical structure 
with outer diameter of about 6 nm, 4 nm pitch, and five 
monomers per turn 0, . The bacterial motility associ- 
ated with type IV pili, called twitching motility, is driven 
by repeated extension, tip attachment, and retraction of 
the pilus filament, by which the bacterium can pull it- 
self forward on surfaces like glass plates or target host 
cells 0. Type IV pili are expressed by a wide range of 
gram-negative bacteria including Myxococcus xanthus 
|29l | and human pathogens Neisseria gonorrhoeae and 
Pseudomonas aeruginosa |30l |. 

The mechanism of retraction is believed to be fila- 
ment disassembly mediated by PilT, a member of the 
AAA family of motor proteins 0, but the microscopic 
details of this process are not known. One might com- 
pare pilus retraction with force generation by micro- 
tubules, which are multistranded filaments with effective 
monomer lengths similar to type IV pili. The helical pitch 
divided by the number of filament strands is 0.8 nm for 
the pili and 0.6 nm for microtubules However, pilus 
retraction generates forces of up to 160 pN H,0|, which 
is an order of magnitude larger than those observed in 
in vitro experiments on microtubules 0, |^ [s^- An- 
other difference is that the pilus retraction velocity is in- 
dependent of filament length Since dissociated pilin 
monomers are stored in the cell membrane waiting to be 
recycled in other filaments j30l |. the implication is that 
the velocity is independent of pilin concentration in the 
membrane. This presumably rules out simple polymer- 
ization ratchet-type models, which have been proposed to 
describe polymerization forces generated by microtubules 
0, 1^ ■ The experimental evidence instead favors a re- 
traction process driven by an active molecular motor Q . 

Pilus retraction is highly processive, and retraction ve- 
locities are of the order of 0.5 — 1 /zm/s 0,|1,E,|33- Gen- 
eration of high forces persists when the PilT concentra- 
tion is reduced, suggestingthat one single PilT complex 
retracts the pilus filament Q . The stall force (the force at 
which the average velocity drops to zero) and the veloc- 



ity at high forces are insensitive to changes in ATP con- 
centration, and the retraction velocity is roughly force- 
independent for small applied forces (< 50 pN) within 
the experimental accuracy 0| . PilT has been shown 
to form a ring structure with sixfold symmetry [s^ , and 
since each subunit has an ATP binding motif, it is pos- 
sible that it can hydrolyze up to six ATP molecules in 
parallel during retraction Q . The outer diameter of the 
ring is about 10 nm, and the inner diameter varies in the 
range 2 — 4 nm [s^ • 

Pilus retraction is interesting from a technological 
point of view, as a potential prototype for a nanoma- 
chine that can generate large forces, and from a biomed- 
ical point of view since pilus retraction is important for 
the infectivity of various severe bacterial pathogens |34l |. 

There are several proposals for how the molecular con- 
stituents of the retraction machinery fit together. One 
of them is that PilT forms a ring around the base of the 
pilus 0, 113 • The hole in the middle of the PilT complex 
seems too small to let the assembled filament through, 
but large enough for pilin monomers. This could allow 
interactions between the pilus and PilT via several ac- 
tive sites (motor subunits) that work together and is the 
principle that we will explore here. For simplicity, we as- 
sume one motor subunit per filament strand and neglect 
possible two-dimensional effects such as angular motion 
of the filament. 

We stress that the purpose of this paper is not to at- 
tempt to describe the detailed molecular mechanisms in- 
volved in pilus retraction, which are largely unknown. 
Rather we examine a new regime of a simple model, 
whose main features arc inspired by experiments. Be- 
low we obtain several results from the model, such as 
large force generation and other properties that agree 
well with interesting experimental results on pilus retrac- 
tion. Moreover, these results are a consequence of corre- 
lations and interactions between the motor subunits and 
are strikingly different from the characteristics of the sin- 
gle building block of the model. 



II. RETRACTION MODEL 

In this section the geometry and equations of motion 
of the model are described. We then discuss the param- 
eters, which come in several kinds: parameters that are 
known for the PilT system, parameters that can be es- 
timated to varying degrees of accuracy, and parameters 
that we will explore in a systematic way. A few param- 
eters cannot be estimated due to the lack of knowledge 
of the molecular details. In this case, we make an ar- 
bitrary choice in order to investigate the qualitative be- 
havior of the model. The basic setup is sketched in Fig. 
^a). A ring of M motor subunits interacts with an M- 
stranded helical filament, with repeat distance d of the 
single strands. The filament coordinate q decreases dur- 
ing retraction. We also allow for deformations of the PilT 
ring. This is described by displacements Xi of the mo- 
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Figure 1: (Color online) (a) The elements of the retraction 
motor model consist of a flexible ring of motor subunits that 
interact with a moving helical filament, (b) Equivalent geom- 
etry after the change of variables yi = Xi — id/M, which places 
the binding potentials of the filament monomers on top of each 
other. The binding potential is assumed to be an asymmetric 
ratchet potential with amplitude U and asymmetry a. The 
undeformed state of the motor protein complex is described 
by the equilibrium positions j/f of the motor subunits, and the 
subunits are elastically confined to their equilibrium positions. 
Due to the helical structure of the filament, the equilibrium 
positions become evenly spread over one period. The mo- 
tor subunits at positions j/i (black circles) interact with the 
filament potential. The open circle represents an unbound 
subunit. During a successful retraction process, a motor sub- 
unit detaches from the filament, relaxes in its confinement 
potential, and rebinds near the next binding site along the 
filament, (c) Distribution of the unbound subunit in B. The 
shaded area represents the probability for the subunit to bind 
to the left of the potential maximum and produce a failed step. 
As the applied force increases, the filament is pulled forward 
relative to the distribution and the probability of a failed step 
increases, (d) Qualitatively new behavior emerges in the limit 
of high stiffness (strong confinement) and many motor sub- 
units. In this limit, all bound motor subunits will not relax to 
the minima in the binding potential. Instead, some subunits 
will interact with the shaded part of the potential and oppose 
force production. In this regime, one expects that the stall 
force depends nonlinearly on the number of motor subunits, 
and nonmonotonically on the stiffness. 



tor subunits from the undeformed state. We assume that 
the motor subunits interact with the filament strands via 
identical one-dimensional binding potentials with period 
d. To mimic the helical structure of the filament, the 
potential of subunit i is displaced a distance id/M rela- 
tive to the potential of subunit 0. The interactions be- 
tween the filament and subunit i therefore have the form 
V{xi — q — id/M), where V is some one-dimensional bind- 
ing potential. We formulate the equations of motion as a 
system of overdamped Langevin equations for q and Xi, 

, , ,dV(xi~q-id/M) I — 

ax, = -kx, - h,{t) — - + ^/2ak^Ut), 



dxi 



M-l 



(1) 



where i = 0, ...,A/ — 1, a and 7 are friction coeffi- 
cients of the motor subunits and the filament respec- 
tively, fc is a spring constant describing the stiffness 
of the PilT ring, hi{t) =0,1 are chemical state vari- 
ables, fee is Boltzmann's constant, T is the temperature, 
F is an external force acting on the filament, and M 
is the number of motor subunits and filament strands. 
Thermal fiuctuations are included through indepen- 
dent Gaussian white noise terms £,i{t) and 'fq(i), which 
obey (6(0) = {Ut)) = 0, = ^^3Kt - t'), and 

{£.q{t)^q{t')) = S{t — t') and have prefactors according to 
the fluctuation-dissipation theorem. The temperature is 
set to T= 310 K. 

Before discussing the terms in the equations of motion 
in detail, we transform the subunit coordinates to place 
the binding potentials on top of each other. The new 
coordinates arc convenient in order to analyze the model 
and are close to previous works with similar models 
I23 . 125I I . The transformed subunit coordinates are yi = 
Xi — id/M + d. An undeformed PilT ring is now described 
as Ui = Hi = d — id/M, which we call the equilibrium 
positions of the subunits. The transformed equations of 
motions are 

ay, = -kiy, - y^) - /i, (i) ^^^-lil + ,/2^ i,{t) 



dyi 



M-l 



dVjy, - q) 
dy. 



(2) 



Note that the binding potentials are all on top of each 
other, as if all motor subimits were interacting with just 
one single strand. The price for this convenience is that 
the equilibrium positions are evenly distributed over one 
period of the potential. The setup in the transformed 
variables is sketched in Fig.mb). 

The motor subunits can be in one of two states: un- 
bound {hi = 0), in which they diffuse around their respec- 
tive equilibrium positions, or bound {hi = 1), in which 
they also interact with the filament through the binding 
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potential V{yi — q). We treat the number of motor sub- 
units Af as a parameter, and we will present results for 
M = 2, 3, 5, 6, and 12. 

For simplicity, we model the binding of the motor to 
the filament by an asymmetric sawtooth potential, shown 
in Fig.^c). For the asymmetry factor we take a = 0.1, 
if not stated otherwise. Some asymmetry is needed to 
give a preferred direction of motion, and the sawtooth 
potential was selected to give a simple parametrization 
of the (unknown) real interaction potential. One possi- 
ble origin of an asymmetric binding potential is surface 
charges on the head of thepilin monomer |3(i| | - for exam- 
ple as described in Ref. Another alternative is that 
the asymmetry can be viewed as an effective description 
of some asymmetry somewhere else in the system, e.g., in 
the direction of the motor steps (power strokes). Based 
on the helical structure of the pilus filament, we take the 
periodicity of the potential to be d = 4 nm 0, |23| , which 
we use for all values of M. 

The amplitude U of the potential is the maximal en- 
ergy required to break the bond between the filament and 
the active site. Pilus retraction is powered by hydroly- 
sis of one or a few ATP per retracted pilin monomer [Sj, 
which sets the energy scale for the potential. Depending 
on conditions, the free energy yield from hydrolysis of one 
ATP in a ceU is about 80 - 100 pN nm ji]. The motor 
subunits are bound together to form the motor complex, 
and we model their confinement with a harmonic restor- 
ing force —kxi = —k(yi —yf). This linear approximation 
is reasonable if the deformations Xi = yi — y^ are small, 
which will be verified below. 

For the binding processes hi{t), we restrict ourselves 
to a sequential reaction scheme with M chemical states. 
We define state j as hk{t) = 1 — 5kj\ i.e., subunit j is 
unbound and the other subunits are bound. The states 
are visited in ascending order, and the (constant) transi- 
tion rate from state j to j + l(modM) is A. The bound 
subunits spend most of their time near a minimum in the 
binding potential, and several geometrical configurations 
are compatible with each chemical state. The main path- 
way for efficient retraction in the model is that the sub- 
units take turns to hop forward to the next minimum as 
they release and rebind to the filament (successful steps). 
When the retraction is less efficient - for example at large 
applied force - the subunits sometimes do not hop for- 
ward (failed step), which leads to geometrical configura- 
tions outside the main pathway. During a successful step, 
a motor subunit goes through the following sequence of 
events: 

(i) The subunit is released from a minimum in the 
binding potential. 

(ii) The released subunit relaxes to its equilibrium posi- 
tion. The filament relaxes in the opposite direction 
due to the forces from the other motor subunits. 

(iii) The subunit rebinds, close to the next minimum in 
the binding potential (otherwise the step fails). At 



the same time, the next motor subunit enters step 
(i). 

(iv) The subunit stays bound and pulls on the filament 
as the other M — I subunits go through steps (i)- 
(iii). After each successful rebinding event, the fil- 
ament retracts a distance d/Ad . 

This mechanism relies on the asymmetry of the poten- 
tial and is similar to the mechanisms studied earlier for 
two elastically coupled particles 13 113 ■ We argue below 
that generation of strong forces in this model relies on a 
binding process that always keeps several motor subunits 
bound to the filament, but the binding order is less im- 
portant. 

We now discuss the parameters of the model. In the 
laser tweezers experiments 0, |^ 0| , the outer filament 
tip binds to an external latex bead with diameter 1 — 2 
fim. Using Stokes law, 7 = QirriR, the approximate 
viscosity 77 = 10^^ pN s/nm^ of the bulk solution sur- 
rounding the cell (somewhere between 10~^ for water and 
8 X 10^*" for glycerin seems reasonable) and R = 1 /im, 
we get 7 « 2 X 10""* pN s/nm for the bead. As a first 
approximation, we neglect the elasticity and friction of 
the pilus filament itself. For the internal friction coeffi- 
cient, we use a = 0.5 x 10"'^ pN s/nm ^ k/X. This sets 
the time scale for internal relaxation a/k much smaller 
than the typical time {M — 1)A~^ between binding and 
release of individual motor subunits and lets the motor 
subunits reach thermal equilibrium between transitions. 
This is consistent with estimates of thermal relaxation 
times over length scales on the order of 10 nm [S^l, which 
is the size of the PilT ring. Another time scale for inter- 
nal relaxation is given by the time to slide down to a 
potential minimum, ad? /U , which we also keep smaller 
than (M — 1)A~"'^. We then expect the velocity to be 
proportional to A, and we will restrict ourselves to this 
quasistatic regime for two reasons. First, this is the bio- 
logically relevant regime where we expect the stall force 
to be independent of A, which corresponds to the experi- 
mental observation that the stall force is independent of 
ATP concentration Second, the exact value of a is 
not critical for the results in this regime, and since it is 
difficult to estimate a accurately, we can avoid making 
our results depend strongly on an unknown parameter. 

Having found useful values for the potential period d 
and the drag coefficients and 7, we go on to investi- 
gate the model behavior as a function of the remaining 
parameters a, M, k, U, and A and properties of the bind- 
ing process. 

III. RESULTS 

A. Methods 

Retraction of the filament means that q decreases, so 
it is natural to study the retraction velocity v = ^dq/dt. 
In the laser tweezers experiment, the tip of the bead is 



held by a static laser trap which is to good approxima- 
tion a harmonic potential - i.e., F — —k^q, with fct on 
the order of 0.1 pN/nm. Numerical solution of Eq. 
using a standard method, known as the Milstein scheme 
[33, produces a deflection trajectory similar to the ex- 
perimental ones. We calculate the retraction velocity by 
fitting a second-order polynomial to a small time interval 
around a point q = —F/kt and take the velocity v(F) as 
the derivative of the polynomial. This is similar in spirit 
to how the experimental data was analyzed H, 0, TTol | . 
The retracted distance —q{t) increases from the initial 
value towards a steady-state, corresponding to the max- 
imal applied force (stall force), which we define as the 
mean applied force in the steady state. To check our 
simulation code, we reproduced analytical results for the 
steady state current in a flashing ratchet model [sl] . 



B. Stall force and force- velocity relation 
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Figure 2: (Color online) Force-velocity relations for high (tri- 
angles), intermediate (squares) and low stiffness (circles) cal- 
culated for M = 5. The y axis is velocity normalized by 
Xd/M, which corresponds to the velocity for M motor sub- 
units if all steps were successful. 



The stall force and force- velocity relation of the motor 
is determined by several competing mechanisms, which 
we now describe qualitatively. 

In the case of a few motor subunits, stalling of the re- 
traction is controlled by two different mechanisms. One 
comes from the finite binding energy between the fila- 
ment and the subunits. Due to the simple shape of the 
potential, we can estimate an upper limit for the stall 
force, using forces AV/Aq. The maximum force that 
each subunit can exert on the filament during retrac- 
tion against an opposing (positive) force is C//(l — a)d, 
so M — 1 bound motor subunits give an upper limit of 
U{M — l)/d{l — a) for the stall force. At finite tempera- 
tures the upper limit is not reached, since the motor sub- 
units diffuse and can pass between potential minima by 
thermal excitation. We think of these thermally assisted 
transitions as slipping events, and they occur more often 
when the subunits are far from their equilibrium positions 
and experience a large confining force. This mechanism 
tends to increase the stall force with increased binding 
strength U. 

The other mechanism has to do with the stiffness k 
and the probability for a step to fail. This probabil- 
ity depends on the distribution of the filament position 
relative to the unbound subunit, which is illustrated in 
Fig. ^c) . The shaded area represents the probability of 
a failed step. Higher probability of failure gives lower 
velocity. Increasing the applied force F at constant stiff- 
ness pulls the filament to the right, increases the frac- 
tion of failed steps, and decreases the velocity. If all 
steps fail, no retraction takes place. The width of the 
position distribution of a subunit in the harmonic con- 
finement is -^/fceT/fc, but the distribution in Fig. Wic) is 
broader, since the filament also fluctuates. The average 
relative position A is roughly proportional to F/k. At 
very low stiffness, the distribution is broad enough for 
some steps to fail without applied load, and the average 
relative position varies strongly with applied force. This 



gives a monotonically decreasing force-velocity relation 
and a low stall force limited by the stiffness. At very 
high stiffness, the distribution is narrow and the average 
relative position is less sensitive to the applied force. Al- 
most no steps fail without applied load, and it takes some 
threshold force before steps start to fail significantly. We 
get a force- velocity relation that is almost independent of 
force at low forces, and the stall force is mainly limited 
by slipping events. Figure |21 shows examples of normal- 
ized force-velocity relations of the model for strong, in- 
termediate, and weak stiffness, compared to the binding 
strength. The curves illustrate the qualitative arguments 
of the preceding paragraphs. We also find that the stall 
force is insensitive to changes in the reaction rate A (not 
shown). Since the reaction rate corresponds to the ATP 
concentration, this is in qualitative agreement with the 
experimental result that the stall force is insensitive to 
the ATP concentration This is expected in the qua- 
sistatic regime that results from the choice of time scales 
discussed at the end of SecE] Stiff systems {k(P/U :s> 1) 
have a plateau in the retraction velocity at low forces, 
which is also the general experimental trend j^, ^3 ■ For 
the parameter regime we have investigated. Fig. |21 gives 
the qualitative shape of the force-velocity relation as a 
fimction of k/U. The size of the plateau is roughly pro- 
portional to the stall force, with a proportionality con- 
stant that is to first approximation a function only of 
k/U. 

As the number of motor subunits and filament strands 
increases, high stiffness can also have a destructive ef- 
fect on the stall force, as illustrated in Fig. ^3. This 
occurs when the distance between equilibrium positions 
is shorter than the region of the binding potential with 
backward slope (shaded) - i.e., M > 1/a. In that case, 
some of the bound motor subunits will tend to interact 
with the shaded region of the potential if the stiffness is 
high. There they act with a negative force on the filament 



and contribute negatively to the force production. This 
effect is enhanced by increased stiffness, and we therefore 
expect the stall force to have a maximum as a function of 
stiffness. This is a qualitatively different behavior than 
with only a few motor subunits (M < f/a) and makes 
the stall force depend in a highly nonlinear way on M . 



C. Parametrization of the stall force 

Equation Q suggests that the stall force might depend 
on the ratio k/U , instead of k and U independently. We 
will use this observation to further analyze force produc- 
tion in the model. A parametrization of the stall force 
is obtained from a combination of the estimated upper 
limit U{M — — a) for the stall force with an func- 

tion of k/U. Using this (unknown) function Jm, which 
also depends on M, we describe the effect of stiffness in 
the following way: 

^^stall(fc, U) = ^^inf /a/(^), Finf = {U - U m) , 

(3) 

where Um is a free M-dependent parameter, independent 
of U and fc, and cP was inserted to make the argument 
of /m dimensionless. One can obtain fuix) by plot- 
ting FstaiiZ-Pinf against k(P /U and adjusting Um- Such 
a plot is shown in Fig. 01 a), with the best fit values of 
Um = 25.4, 25.9, 28.9, 28.3, and 33.5 pN nm for A/= 
2, 3, 5, 6, and 12, respectively. For A/=3, 5, and 6, the 
data points fall on a single curve for each M to good ap- 
proximation, while M—2 and 12 show some scattering. 
We can understand the ansatz Q and the results in Fig. 
13 in light of the qualitative arguments in Sec. IIII 51 We 
interpret i^inf as the stall force in the stiff limit fc ^ oo 
and Um as the effect of thermal fluctuations, which in- 
duce slipping events and thereby lower the maximal force 
that the binding potential can support. The stiff limit 
k oo can be simulated by locking the motor subunits 
to their equilibrium positions. The results arc shown in 
Fig. |3Jb), using the above values for Um- The stall force 
per bound motor subunit is consistent with Eq. Q and 
liiiix -t fM{x) = 1 for M = 2, 3, 5, and 6. The case 
M = 12 also gives a straight line for the stall force as 
a function of U in the stiff limit, but neither consistent 
with lim^j^oo fi2{x) = 1 nor with the value of U12 from 
finite stiffness. This illustrates the effect of dense motor 
subunits {M > 1/a) in a stiff system. 

The function fM{kd?/U) is a normalized stall force per 
motor subunit. As is evident from Fig. |3fa), the force 
production per motor subunit varies strongly with stiff- 
ness and also with the number of subunits, M . Clearly, 
the interaction between the motor subunits is impor- 
tant for the force generation. The stall force is more 
complicated than a sum of contributions from the indi- 
vidual parts. According to the arguments in Sec. IIII Bl 
one expects /m^x) to be an increasing function of x for 
M < 1/a, since larger stiffness decreases the probabil- 
ity that the unbound motor subunit binds to the wrong 
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Figure 3: (Color online) (a) Stall force as a function of k/U 
in the range 5 < fc < 250 pN/nm and 80 < U < 200 pN nm. 
The normalization Fini is chosen according to Eq. Lines 
are guides to the eye. The case of many subunits (Af — 12) 
is qualitatively different from the other cases. Here, the stall 
force depends nonmonotonically on k/U. (b) Stall force per 
bound motor subunit in the stiff limit fc — > 00 as a function 
of L'^ — Um. For Af < 6 the simulations are consistent with 
/a/(cx)) = 1 in Eq. @ (solid line). M = 12 gives a lower 
stall force that does not fit Eq. which again illustrates 
the nonlinear behavior of the stall force on AI. Error bars in 
both graphs are smaller than the symbols. 



potential well. For M > 1/a, the stall force is expected 
to have a maximum as a function of stiffness, reflecting 
the inset of the qualitatively different behavior at high 
stiffness. These expectations are consistent with Fig. O 
where f2{x), faix), /a (a;), and fe{x) are monotonically 
increasing in the simulated region, whereas fi2{x) has a 
maximum around 1.5. To confirm that dense motor sub- 
units give lower stall force also for lower values of M, we 
performed simulations with a = 0.3 and found that with 
U = 200 pN nm and Af = 5, a maximum stall force oc- 
curs for stiffness between 40 and 150 pN/nm (not shown), 
as expected. 
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F[pN] 

Figure 4; (Color online) Force-velocity characteristics from 
simulation and experiment. Triangles; average experimental 
velocity from Fig. 3(b) of Ref. Thick solid curve: A4 = 5 
model with U = 200 pN nm, fc = 25 pN/nm, F^taii = 123 pN, 
and A = 1180 s~^. Diamonds: a single retraction event from 
Ref. 13. Dashed curve: M = 5 model with U = 200 pN nm, 
= 60 pN/nm, F^tM = 160 pN, and A = 1800 s^\ Isolated 
motor subunits - i.e., simple flashing ratchet models (black 
symbols) are qualitatively different, as illustrated for M = 1, 
U = 200 pN nm, A = 2000 s~\ and toff = 25 ms (*) or 2.5 
ms (*). 



D. Comparison with an individual motor subunit 



To further highlight the effect of interactions and cor- 
relations, it is interesting to compare the interacting ring 
model to an isolated motor subunit. This is a single par- 
ticle in a flashing ratchet potential and a special case of 
the interacting model, with M = 1 and k — > oo. This 
single-particle flashing ratchet (SPR) has been studied 
extensively 0, [s^ in different versions, and its properties 
are qualitatively different from those of the interacting 
model in several interesting respects. For an SPR in this 
simple version, the mechanism to pull the particle to the 
next potential well is not present. Instead, forward mo- 
tion relics on thermal noise to make the particle diffuse 
forward while in the unbound state. This means that 
the probability for a forward step can never exceed 1/2 
for a single chemical cycle even without applied force. 
With nonzero applied force, the free diffusion is superim- 
posed on a backward motion with velocity Wdrift = —F/l] 
hence, the velocity is substantially reduced even at very 
low forces. These features conspire to make both the 
maximal velocity and the stall force depend strongly on 
the friction constant, the reaction rate, and how much 
time the particle spends unbound during a reaction cy- 
cle [l^^ • Two examples of force- velocity relations for 
SPR are shown in Fig. 21 with 4 nm periodicity, U = 160 
pN nm, A = 2000 s~^, and damping as above. The tran- 
sitions between the bound and unbound states are deter- 
ministic in these cases, and the time toff spent unbound 
during each cycle was 2.5 ms and 25 ms. 



E. Application to a real motor system 

In Fig. 21 we compare model results with experimental 
data on N. gonorrhoeae. Several force-velocity relations 
averaged over many retraction events, as well as two sin- 
gle events, are preesented in Ref. H,[3. The maximum 
forces produced in the different events are distributed 
between 40 and 160 pN. There are several experimental 
factors that can make the maximum measured force in a 
particular retraction event lower than the intrinsic stall 
force 0,jlOj|. These factors include, for example, breakage 
of the filament and are not included in the model. Ac- 
cordingly, the model should be compared to the data that 
reach the highest forces. This leaves one single trajectory 
0, shown in Fig. ^ and several average curves. The av- 
erage force- velocity data in Fi g. |^ represent an average 
of the data in Fig. 3 of Ref. [ifl- Both experimental 
curves give positive velocities up to about 160 pN, which 
is in the upper tail of the maximum force distribution 
[To| . but the velocity at lower forces differs considerably 
between the two curves. 

To describe the PilT system, there are two natural 
choices for M: namely, M = 5, reflecting the five-fold 
symmetry of the filament, and M = 6, reflecting the 
six-fold symmetry of the PilT molecule. We present sim- 
ulation results for M = 5, but 7\/ = 6 is qualitatively no 
different. 

The model can describe both experimental curves to 
some extent. The parameter values that describe the 
single event and the average data (see Fig. 2)) differ in 
reaction A rate and stiffness k. There are two important 
deviations, which we discuss next. 

The average velocity falls off exponentially at high 
forces while the results of the model decay faster. 
The exponential decay can be described by an Arrhcnius 
law for the rate limiting step The present model 

does not account for this behavior. 

The single retraction event has a different decay at 
high forces and agrees better with the characteristics of 
the model. However, the single event shows an initial 
increase in velocity at low forces, rather than a plateau. 
Such initial increases in velocity are also obtained in some 
individual simulation runs, but disappear when the av- 
erage is computed. All the experimental data also suffer 
from a possible systematic underestimation of the veloc- 
ity near F=0 H,0|- Within the experimental accuracy, 
the general trend is a constant velocity up to about 50 pN 
Hi EOl I which is consistent with the results of the model. 

F. Deformations and elastic approximation 

From the results presented in Figs. 01 anc^l and Eq. 
©, it is possible to estimate the magnitude of the de- 
formations Ui of the motor complex near the stall force. 
This is useful, since the harmonic approximation for the 
confining force on the motor subunits is questionable for 
large displacements. At the stall force, with il/ — 1 sub- 
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units bound to the filament, the mean displacement of 
a subunit can be estimated to Fsta\\/k{M — 1). As seen 
earlier, the stall force is less than (A/ — l)U/d{l — a), 
which gives a displacement less than U /kd{l — a). From 
Fig. 121 the interesting and relevant regime with high nor- 
malized stall force has kd/U > 1 nm^^. This gives a 
displacement less than 1 nm near the stall force. This is 
not excessively large compared to the dimensions of the 
PilT ring, whose diameter is 11.5 nm [3^ . 

G. Role of order in the binding process 

We next examine the effect of the ad hoc assumption of 
a sequential binding process on the results of the model. 
We compared the ordered sequence with two less corre- 
lated binding schemes. For random order with only one 
subunit free at the same time, the stall force is essen- 
tially unchanged, but the velocity decreases with about 
50%. Alternatively, if the binding and unbinding events 
are assumed to be independent for the different subunits, 
all subunits might occasionally become unbound simul- 
taneously, which releases the filament from the motor. 
The resulting retraction events become highly irregular 
and have low stall force and mean velocity. We conclude 
that some degree of correlation between the motor sub- 
units is necessary in order for the model to simulate the 
experiments. 



IV. DISCUSSION AND CONCLUSION 

We investigate a model for force generation in finite en- 
sembles of motor subunits interacting through an elastic 
backbone, which is inspired by the pilus retraction ma- 
chinery, the strongest molecular motor reported so far. 
The model is prototypical, rather than realistic in detail, 
and offers a possible mechanism for generation of large 
retraction forces. It includes a ring of motor subunits 
surrounding the pilus filament, following a suggestion in 
Ref. |35l |. Some parameters in the model can be roughly 
estimated, and we explore parts of the remaining param- 
eter space and focus on generation of large forces. 

We find that the stall force depends on the binding 
strength U between motor and filament, the stiffness k 
of the motor complex, the number of motor subunits M, 
and an asymmetry parameter a. For high enough stiff- 
ness we find qualitatively different properties compared 
to the well-studied model of a single particle in a flash- 
ing ratchet potential, which is the basic building block of 
our model. This is not surprising, since the mechanisms 
that generate motion are different in the two cases. The 
motion in the flashing ratchet model is dependent on dif- 
fusion Q. Our model also contains diffusive motion of 
the motor subunits, but diffusion is not necessary for the 
motor to work j2J, |23 . The dependence of the stall force 
on U, k, and Af is well parametrized in empirical scaling 
plots (Fig.O. The scaling ansatz in Eq. Q relics on the 



presence of an interaction between the motor subunits, 
and the stall force depends nonlinearly on the number of 
motor subunits. 

Low stiffness compared to binding strength has a 
strong destructive effect on the force production. For a 
small number of motor subunits, the stall force increases 
monotonically with increasing stiffness, but when the mo- 
tor subunits become dense enough there is a crossover to 
a different regime, where a high stiffness instead has a 
destructive effect on the stall force. For the asymme- 
try studied here, the crossover occurs around M = I /a. 
For M = 5 or 6, corresponding to the number of flla- 
mcnt strands or PilT subunits, a system with a > 0.2 
would have the interesting property that the maximum 
stall force is obtained for finite stiffness. If the binding 
strength between the PilT subunits can be genetically 
engineered, it might be possible to observe this effect in 
future experiments. 

We compare results from the motor model for the fila- 
ment retraction force- velocity characteristics with exper- 
iments on N. gonorrhoeae |a, llo| . Since the molecular 
details of the retraction mechanism are unknown, it is 
unclear to what extent the agreement we see refiects ac- 
tual similarities between the model and the real system. 
More information would be useful for the construction of 
more detailed models. Nevertheless, the model can de- 
scribe the general features of the experimental results - 
i.e., the plateau in the force- velocity relation at low ap- 
plied forces and the high stall force that is independent 
of reaction rate 0, ^3 ■ 

For a quantitative comparison we select two differ- 
ent experimental force-velocity relations that show large 
forces, one single event, and one averaged curve, as shown 
in Fig. 21 With different parameters, the model gives 
a reasonable description of the average data, as well as 
what looks like an atypical single event. The model devi- 
ates from the average data at high forces, which indicates 
that something is missing from the description. However, 
the data might include variations in cellular conditions 
that affect the average at high forces. This is not ac- 
counted for in the simulations, where the average is taken 
over thermal fiuctuations with fixed parameters. The sin- 
gle event is described better and does not suffer from such 
a complication. In this case, the model is limited at high 
forces mainly by thermally assisted transitions of motor 
subunits between potential minima. Given the simplic- 
ity of the model, we find the agreement with experiments 
encouaging. 

Small ensembles of interaction motor systems have pre- 
viously been found to possess rich behavior without ap- 
plied force p3 . l25ll27l |. As we have shown, this is true also 
in the limit of high forces, which is a realistic experimen- 
tal situation. Generation of strong forces in nanoscale 
devices is also of technological interest, and it is tempt- 
ing to speculate about the possibility to realize a setup of 
interacting motors units pulling on an artificial filament 
such as a carbon nano tube. 
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